The negative index of refraction demystified 



Martin W McCallf|, Akhlesh LakhtakiaJ and 
Werner S Weiglhofer§ 

f Department of Physics, The Blackett Laboratory, Imperial College of Science, 
Technology and Medicine, Prince Consort Road, London SW7 2BZ, Great Britain 
| CATMAS — Computational & Theoretical Materials Sciences Group, Department 
of Engineering Science and Mechanics, Pennsylvania State University, University 
Park, PA 16802-6812, USA 

§ Department of Mathematics, University of Glasgow, University Gardens, Glasgow 
G12 8QW, Great Britain 

Abstract. We study electromagnetic wave propagation in mediums in which the 
effective relative permittivity and the effective relative permeability are allowed to 
take any value in the upper half of the complex plane. A general condition is derived 
for the phase velocity to be oppositely directed to the power flow. That extends the 
recently studied case of propagation in mediums for which the relative permittivity 
and relative permeability are both simultaneously negative, to include dissipation as 
well. An illustrative case study demonstrates that in general the spectrum divides into 
five distinct regions. 
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1. Introduction 

Materials with either negative real relative permittivity or negative real relative 
permeability at a certain frequency are not uncommon. Consideration of materials for 
which both these quantities are simultaneously negative and real-valued, commenced 
with Veselago's paper of 1968 Although he pointed out many unusual properties 
of such materials, including inverse refraction, negative radiation pressure, inverse 
Doppler effect, over three decades ago, the considerations were completely speculative 
in view of the lack of a material, or even a nonhomogeneous composite medium, with 
a relative permittivity having a negative real part and a very small imaginary part. A 
breakthrough was achieved by Smith et al. ||, who, developing some earlier ideas by 
Pendry et al. 0, £|, ||, presented evidence for a weakly dissipative composite medium 
displaying negative values for the real parts of its effective permittivity and effective 
permeability. Their so-called meto-material consists of various inclusions of conducting 
rings and wires embedded within printed circuit boards. Their conclusions were based 
on observations from three separate composite mediums. 

• Medium 1 consisted of a lattice of ring-like inclusions, which for a certain field 
configuration was presumed to have a resonant relative permeability fi e ff(u) of the 



where the resonance frequency ooq depends principally on the geometry of the rings. 
In this model, dissipation is facilitated by T, and F (0 < F < 1) is the ratio of the 
area occupied by a ring and that of a unit cell. For weak dissipation, the real part 
of [j, e ff is negative for u < uj < u / y/1 — F. 

• Medium 2 consisted of an included matrix of wires. The effective relative 
permittivity e e ff of this composite medium supposedly displays plasma-like 
behaviour according to 



the effective plasma frequency u p being dependent on the geometry. In such a 
medium, e e ff is negative for uj < u> p . 

• Medium 3 combined the first two, it being postulated that the combination would 
exhibit negative real-valued permeability and negative real-valued permittivity 
within a certain frequency range. Although some numerical and experimental 
evidence was presented in 2000 , the most striking report appeared the following 
year || and gave preliminary indications of the observation of the so-called negative 
index of refraction. 
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Other types of nanostructural combinations with similar response properties can also 
be devised 0]. 
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The emergence of a clear terminology is often a difficult process with regards 
to scientific findings relating to novel effects, something that is also apparent in the 
present instance. The original classification of the materials exhibiting the effects 
discussed labelled them left-handed electromagnetic materials But chiral materials 
are important subjects of electromagnetics research and the terms left-handedness and 
right-handedness have been applied to the molecular structure of such materials for 
well over a century ||. The continued use of the term left-handed materials (LHMs) 
for achiral materials in, for example, fj, ^, f| will thus confuse the crucial issues. Other 



authors |10| are promoting the term backward medium (BW) which presumes the a 
priori definitions of forward and backward directions. In the most recent contribution 
|TT| (which also provides the most extensive theoretical and numerical analysis of 
the negative index of refraction to date), the authors introduce the technical term 
double negative (DNG) medium to indicate that the real parts of both permittivity 
and permeability are negative. While sensible enough, such nomenclature conceals the 
importance of dissipative effects. 

In time, a consensus about terminology will undoubtedly emerge; and it is not 
our aim to contribute to this particular discussion. Instead the purpose of this note is 
pedagogical. In the first instance, it is important that dissipation be included in the 
analysis. This was largely neglected in the reports cited earlier, with the exemption of 
the most recent study fllf . 

Secondly, it is desirable to derive the general condition for the type of anomalous 
propagation that is characteristic of the considered materials: namely, where the phase 
velocity is directed oppositely to the power flow. 



2. Plane Wave Propagation 



Consider a plane wave propagating along the z axis in a linear, homogeneous, isotropic, 
dielectric-magnetic medium whose relative permittivity and relative permeability are 
denoted by e r and /i r , respectively. An exp(— iut) time-dependence is assumed here. 
Then 



E(z) = A exp(ik nz) u x 



nh 



Biz) = — V x E(z 

icu uj 

Ti 

H(z) = A exp(ikonz) u 



A exp(ikonz) u 



y ■ 



v ■ 



(3) 
(4) 
(5) 



where k is the free-space wavenumber, r] is the intrinsic impedance of free space, and 
n 2 = e r fi r . Consequently, the Poynting vector is parallel to the z axis and its time- 



average is given as 
P,(n) 



I 

— u, 

2 



Re 



E(z) x W(z) 



Re 



n 



\A\ 2 
2^ 



exp(— 2fcoIm[n]z) , (6) 



where Re[-] and Im[-], respectively, denote the operations of taking the real and the 
imaginary part, whilst * indicates complex conjugation. 
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Let us now assume a Lorentzian model for e r and /i r . This will include the specified 
forms (0) and (0) as special cases. Dissipation results from the imaginary parts of e r 
and \i r whilst causality dictates that Im[/i r ] > and Im[e r ] > 0, so that e r and \i r lie in 
the upper half of the complex plane. 

However, there are two resultant complex refractive indexes, n± = ±^e r fi r , of 
which n + lies in the upper half of the complex plane and n_ in the lower half. The 
situation is summarized in Figure 1 for which the resonant form of (H) was used as 
representative of both e r (u) and fJ> r (uj). Of course, the resonances of e r (ui) and fJL r (u) 
are unlikely to coincide, so that for a particular value of u, the arrows corresponding to 
e r and \i r will not be necessarily parallel. Only the upper half of the complex plane is 
shown in the figure. 

Now n± may be written as 

n± = ±n exp i<p n , (7) 

where 



n 



|/i 



+ 



(8) 



Here 4> e and representing the arguments of e r and pi r respectively, must obey the 
conditions < 4> €jfl < ir. Consequently, < n < tc. We then always have 



and also 



Rc 



Re 



I'r 



71- 



> i.e. P z {n + ) > 



<0 i.e. PJn)< . 



(9) 



(10) 



Thus the choice n + always relates to power flow in the +z direction, whilst n_ always 
relates to power flow in the — z direction. Since necessarily Im[n + ] > and Im[n_] < 0, 
power flow is always in the direction of exponential decrease of the fields' amplitudes. 

We can now identify when the phase velocity is opposite to the direction of power 
flow. This occurs whenever Re[n + ] < (and consequently Re[n_] > 0, also). After 
setting 



e' + ie" 



(where e' r ,e" and fi' r , fi" are the real and imaginary parts of the relative permittivity 
and the relative permeability, respectively), the following condition is straightforwardly 
derived for such propagation: 



+ {<■ 



//2\l/2 



- e r 



i / / 2 //S 



1/2 



(12) 



Before turning to a fully illustrative example in the proceeding section, let us 
investigate some immediate repercussions of the inequality fll2"|) which is central to this 
paper. 
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Consider, in the first instance, the behaviour at a resonance of the relative 
permittivity, i.e. e' r = 0, e" > 0. Then, ([12]) reduces to 

</2\l/2 



>/4', (13) 



an inequality that is always fulfilled when fj,' r < 0. 

Likewise, at a resonance of the relative permeability, i.e. fi' r = 0, fi" > 0, fljjD is 
fulfilled whenever e' r < 0. 

Further insight into inequality (0) can be gained by requiring that 

+ (e' r 2 + ef) 1/2 >e' r + e>: (14) 

and 

. / /2 . //2\V2 , „ 

+ I /i r + /i r J >jj, r + /j, r (15) 



simultaneously hold. Consequently, ( |I2]) is definitely satisfied. It should be 



remarked though, that the parameter space of the permittivity and permeability 
that fulfils ( |14D and (|T5| ) is only a subset of the one fulfilling (|T2"D. In any case, 



( |T4"D holds if and only if e' r < 0, and ( |T5"D holds if and only if < (we remind 
the reader that e" > 0, /jl" > because of causality requirements). We note that 
e' r < and fi' r < can only occur close to absorption resonances (as discussed in 
the previous item). 

Finally, consider an electromagnetic wave propagating in a plasma below the plasma 
frequency (e' r < 0,fi' r = 1) and in which dissipation is very small (e'f 1, /jl" <ti 1). 
Straighforward Taylor expansions reduce inequality ([12]) to 



e" 



141 >d- as) 

/'r 

Therefore, the existence of the type of anomalous propagation being studied here 
depends in this case crucially on the ratio of the imaginary parts of the relative 
permittivity and relative permeability. Whether the criterion (|i~6|) is satisfied or 
not, the power flow in this case is in any case small. 

3. A detailed illustrative case study 

Let us exemplify the foregoing in detail by an explicit invocation of the Lorentz model 
for both e r and /i r ; thus, 

*•(*<>) = 1 + ~, ~ ^2 > (17) 

1 + [N- 1 - zA e A 1 



/V(A ) = H ; —2 



-1 



1 + [K-n 1 - ^A m A " 

Here Ao = 27r/fc is the free-space wavelength, p £tm are the oscillator strengths, 
A e ,m(l + ^(Tm) -1 ^ 2 are the resonance wavelengths, while \ e ,m/N e>rn are the resonance 
linewidths. 



The negative index of refraction demystified 



6 



Figures 2a and 2b comprise plots of the real and imaginary parts of e r and \i r 
as functions of A , when p e — 1, p m — 0.8, N e = N m = 100, A e = 0.3 mm and 
A m — 0.32 mm. Clearly, five separate spectral regions can be identified in Figure 2. 
At the either extremity of the horizontal axis are the two regions wherein e' r > and 
/i', > 0. In the neighbourhood of Ao = 0.22 mm, e' r < but fi' r > 0. Both e' r < and 
/i', < in the neighbourhood of A = 0.25 mm. Finally, e', > but fi' r < around 
A = 0.31 mm. Of course, both e" > and [i!' r > for all A . 

Detailed calculations confirm that the spectral region wherein the inequality ([12]) is 
satisfied is larger than the middle region (wherein both e' r < and /i', < 0). The former 
cover parts of the adjoining regions, in which either (|14]) or (|l|) holds. 

In the five spectral regions identified, the isotropic dielectric-magnetic medium 
would respond differently to monochromatic electromagnetic excitation. Suppose that 
a plane wave is normally incident on a half-space occupied by this medium. The 
reflectance -R(A ) is then given by the standard expression 



where < R < 1 for all Ao by virtue of the principle of conservation of energy. The 
reflectance spectrum calculated with the constitutive parameters used for Figures 2a and 
2b is shown in Figure 3. The reflectance is markedly high in the two regions wherein 
e' r and fjf r have opposite signs, but not in the other three regions. The reflectance is 
particularly low in the leftmost and the rightmost regions (e' r > and fi' r > 0) because 
the ratio fi r (X ) / e r (\o) is close to unity therein. However, the reflectance is somewhat 
higher in the central region (e', < and [i! r < 0) because |/x r (A ) / e r (A ) | < 0.25. 

4. Conclusions 

In this pedagogical note, we have derived a general condition for the phase velocity to be 
oppositely directed to the power flow in isotropic dielectric-magnetic mediums in which 
the only constraints on the values of the relative permittivity and relative permeability 
are those imposed by causality. In this regard, the topical case of mediums in which e r 
and \i r have negative real parts is seen to be a sufficient, but not necessary, condition 
for such propagation, as noted in the comments succeeding (|12|). An illustrative case 
study has shown that there are, in general, five distinct spectral regions, characterized 
by the various sign combinations of the real parts of e r and [i r . 
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Figure Captions 



Figure 1. Argand diagram parametrically displaying fi e ff(cu) from equation ([!]), with 
T = 0.1u and F = 0.5. On taking (|l|) as a model resonance form for meta-materials, 
the plot can also be regarded as displaying the effective permittivity e e //, for which 
the resonance at uj = ujq is unlikely to coincide, and hence the arrows indicating e e ff 
and fi e ff do not coincide in general. The complex number n + = -\--J^ e ff^effi while the 
corresponding index n_ = —y/^^JJJLTff is not shown. The dots indicate equi-spaced 
frequencies from u = to u = 2ujq. 

Figure 2. (a) Real parts of the relative permittivity and relative permeability according 
to equations fllTD and (]T8|) , respectively, when p e = 1, p m = 0.8, N e = N m = 100, 
A e = 0.3 mm and A m = 0.32 mm. The significance of the identified five regions of the 
spectrum is explained in the text, (b) Imaginary parts of the relative permittivity and 
relative permeability according to equations (0) and fll8|) , respectively, when p e = 1, 
p m = 0.8, iV e = N m = 100, A e = 0.3 mm and X m = 0.32 mm. 

Figure 3. Plane wave reflectance -R(Ao) calculated with the constitutive parameters 
depicted in Figures 2. 
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